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Abstract 

Exact solution of Dirac equation for a particle whose potential energy and mass 
are inversely proportional to the distance from the force centre has been found. The 
bound states exist provided the length scale a which appears in the expression for the 
mass is smaller than the classical electron radius e 2 /mc 2 . Furthermore, bound states 
also exist for negative values of a even in the absence of the Coulomb interaction. 
Quasirelativistic expansion of the energy has been carried out, and a modified expres- 
sion for the fine structure of energy levels has been obtained. The problem of kinetic 
energy operator in the Schrodinger equation is discussed for the case of position- 
dependent mass. In particular, we have found that for highly excited states the 
mutual ordering of the inverse mass and momentum operator in the non-relativistic 
theory is not important. 
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1 Introduction 



The concept of the effective mass in theoretical physics is quite efficient because it allows 
to reduce a many-body problem to a single particle one, without the loss of the main 
contributions into the mechanism of the formation of various physical phenomena from the 
inter-particle interactions. We can exemplify this statement by specific problems from the 
semiconductor physics, superfluid 4 He theories, problems of nanophysics as well as by a 
number of other problems from the condensed matter physics PjJ El El HI El ■ At the same 
time, the suggested approach raises the issue of mutual ordering of momentum operator 
in the Schrodinger equation and the inverse effective mass in the kinetic energy operator, 
which is the momentum function of the particle coordinate El HE GUI ^2 • This problem 
seems to disappear if one proceeds from the Dirac equation, but the transition to the 
non-relativistic case is far from being simple as it might look at the first glance 

In this paper we suggest a solution of the Kepler problem (i.e., a study of the particle's 
movement in the Coulomb potential) in the Dirac theory for a particle with the given 
effective mass m* dependent on coordinate r. We assume this dependence to be relevant at 
distances of the particle from the force centre smaller than the Compton length A = h/mc, 
where m is the mass m* at r — > oo. For such distances the notion of the particle coordinate 
is lost as attempts to localize the particle in the space domain with linear dimensions ~ A 
lead to the creation of new particles. In this connection, it seems possible to try and take 
effectively into account the processes of the interaction of the particles with the force centre 
at the subatomic scales through a coordinate-dependent mass. 

In the general case, we assume that, e.g., for an electron the value of m* can be presented 
in the form of a multipole expansion 



where a, b, . . . are constants formed by the mechanism of the particle interaction with the 
vacuum fluctuations in the presence of the force centre; later on we consider these constants 
to be the given initial parameters of the problem. The processes of the interaction at the 
subatomic scales within the quantum field theory lead, in particular, to the deformation 
of the Coulomb interaction in the atom when the distance between the electron and the 
nucleus is small. In other words, it is quite possible to consider the electron charge as a 
function of the particle coordinate. We can therefore speak about the effective considera- 
tion, in the Dirac theory, of the radiation effects which are due to the renormalization of 
the electron mass and charge; that is why we can make an attempt to account, within such 
a phenomenological approach, for the observable superfine structure of the atom energy 
spectrum. 

The above-mentioned arguments can be disregarded. Then, one treats the Dirac equa- 
tion formally as a mathematical problem in which the mass m* is dependent on the radius- 
vector r. 

Below, we provide an exact solution of the Kepler problem in the Dirac theory for the 
case when only the first two terms are taken into account in the expansion for the effective 
mass (jl.Ij) . 
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2 Initial equations 

We start from the Dirac equation (in familiar notation), i.e., 



[a.f>) c + m*c 2 (3 + U 



tp = Eif), 



where d, (3 are the Dirac matrices, the particle effective mass is 

* ( i a ' 
m — m 1 + 



(2.1; 



(2.2) 



and the energy of the Coulomb interaction is 



U = , 

r 



(2.3) 



where e is the particle charge. 

In the following we calculate the energy spectrum E. To do this, we introduce the 
function ijj so that 



ap) c + m*c 2 (3 + (E - U) 



(2.4) 



that is, we apply the substitution used also in [T2"l ITB*] . For the new function ip from the 
Dirac equation (|2.1j) we find 



ap) c + m*c 2 (3 - c [(dp) U-U (dp)] - (E - U) 2 } ip = 0. 



It is easy to see that the commutator 



[dp, U] = —ih(crn)/3 



,dU 
dr ' 



and 



dm* 



[(ap)c + m*c 2 /3] 2 = p 2 c 2 + m* 2 c A -ihc 3 ((rn){3" , . 



dr 



where <r = (a x , a y ,cr z ) are the Pauli matrices, 



-/ 

/ o 



are 4x4 matrices, n = r/r is a unit vector. 

Our initial equation ()2.5|) now takes the following form 



p 2 c 2 + m *2 c 4 _ ihc 3 {(Tll) p» ^! + iTic (tm)p! ^-(E- Uf U = 0. 

dr dr 



(2.5) 



(2.6) 



(2.7) 



(2.8) 



Using explicit expresions for m* and U (see Eqs. (|2.2j) (|2.3Jl ). we obtain, after simple 
calculations, 



2m 



E 



mc 2 a 



1 1 

r 2mr 2 



ih 



— (an){mc 2 a(3" + e 2 /?') + m 2 cV - — 



£ 2 ™ 2 - 4 



m c 



2mc 2 



(2.9) 



We note that our equation has the form of the Schrodinger equation for a particle moving 
in the Coulomb potential with an addant to the centrifugal energy. 



3 The radial equation 

We now pass to spherical coordinates in Eq. (|2.9|) : 



h 2 1 d 2 1 



r + 



2m r dr 2 2mr 
Et 2 2 \ 1 1 £ 2 - /,rc ! - 



f 2 , ^/ W 2 a" i 2 A/\ , 2 2 2 e 

L H (crnjfmc ap +ep)+mca r 

c c 2 



mc 2 / r 



m ° a I - ^ = 2mc 2 ^' (3,1 ' 



where L is the angular momentum operator. 

As the operator in square brackets depends solely on the angles, the variables in Eq. (jrS.l}) 
can be separated. We further make use of the fact that 

L 2 = (trt)[(trL) + h], 

and introduce the following operator 

A = -[(<rL) + h] + Z -((rn)(mc 2 ai3" + e 2 /3'), (3.2) 

with 

A 2 = [(<TL) + h} 2 + m 2 c 2 a 2 -e 4 /c 2 . 

It is now not difficult to show that the operator in square brackets in Eq. (|3.1j) equals 
A(A + h). Then, we rewrite Eq. (13. 1|) as follows 

h 2 1 d 2 A(A + h) ( E 2 2 V \ 7 E " ~ mV 7 ^ 

~ o + ~^ 2 — 2 e -mca)~)ijj= — - — - — ip. 3.3 

2mrdr z 2mr z \mc A J r J 2mc z 

As (<tL) = (J 2 — L 2 — S 2 )/ft, where the total angular momentum J = L + S, and the 
spin operator S = ha/2, it is easy to see that the operator 



j2_ L 2_g2 \ 2 / e 2N2 



A 2 = | — + h\ + (mca) 2 



It then follows that the eigenvalues of this operator 

A 2 = h 2 



2 / 2\ 2 

+ (mca) 

c 



j = l±l/2, 1 = 0,1,2,.... 
Further, if j = I + 1/2, then 

A 2 = h 2 (l + l) 2 + (mca) 2 - (e 2 /c) 2 , 

and if j = I — 1/2, then 

A 2 = h 2 l 2 + (mca) 2 - (e 2 /c) 2 . 
Now it is not difficult to find the eigenvalues of the operator A ()3.2|) : 



A = -hy/(l + l) 2 + (mca/h) 2 - (e 2 /hc) 2 , 

for j = / + 1 /2 and 

A = hy/l 2 + (mca/h) 2 - (e 2 /hc) 2 , 
for j = I — 1/2. Finally, we can easily find the eigenvalues of the operator A(A + K): 

k(k + h) = h 2 i*(i* + i), 

where the quantum number 



/* = y/(j + 1/2) 2 + (mca/h) 2 - (e 2 /hc) 2 - 1/2 t 1/2. (3.4) 

In Eq. (13. 4 j) . the upper sign is for j = I + 1/2 while the lower sign is for j = I — 1/2. 

Now, substituting the operator A (A + K) with its eigenvalues in Eq. (|3.3j) . we can write 
down the equation for the radial part R of 

/i 2 1 d 2 fr 2 Z*(r + l) f E 2 2 \ll D ^ 2 -m 2 c 4 D 

r^ r H — 5 — - - — ^ e ~ mc a I -f R = o — R - ( 3 - 5 ) 

2mrdr 2 2mr 2 \mc 2 J r ) 2mc 2 

4 Energy eigenvalues. Discussion of the results 

Formally, equation ()3.5j) coincides with the non-relativistic Schrodinger equation for the 
Kepler problem with the energy 

E 2 — m 2 c 4 
2mc z 

and the charge squared 

mc 2 



For the existence of the bound states, it is necessary that the "potential energy" in 
Eq. ()3.5|) should be of the attractive nature or, otherwise said, the charge squared e* 2 
should be a positive value e* 2 > 0. 

Then one can write the Bohr formula for the energy levels E*\ 



E* 



me 



* i 



2fr 2 (n r + Z* + l) 2 ' 

where n r — 0, 1, 2, . . . is the radial quantum number, hence, the equation for E 

E 2 — m 2 c 4 m(Ee 2 /mc 2 — mc 2 a) 2 
2mc 2 2h 2 n* 2 

where we introduced the "principal" quantum number 

n r + I* + 1 = n r + + 1/2) 2 + (a/A) 2 - a 2 + 1/2 t 1/2, 



n 



where a = e 2 /hc is the fine structure constant and A = h/mc is the Compton length. 
Solving this equation we finally find for the energy spectrum 



E 



mc 



1 + (a/n*y 



a a 
A n* 2 



a" 



(a/A)' 



n 



*2 



(4.2) 



The wave functions R are the usual radial functions of the non-relativistic hydrogen 
problem with the quantum numbers n*, I* and with the charge squared e* 2 . In order to 
determine the full wave function ip it is necessary to substitute the function ip into Eq. ()2.4j) 
in the form of the product of the radial function R and the spherical spinor. 

The bound states exist when e* 2 > 0. Using Eq. (j4.1j) for the energy E in Eq. ()4.2|) we 
find the condition for the parameter a that does not depend on the quantum number n*: 
a < e 2 /mc 2 . Thus the bound states exist for arbitrary negative values of the parameter a. 
If the length parameter a > 0, it should be smaller than the electron classical radius. In 
other words, these are the distances 1/a ^ 137-fold smaller than the Compton length A, 
where the well-defined notion of the particle coordinate is lost. It is interesting to note that 
for a = e 2 /mc 2 the energy E = mc 2 , hence, there is just one level. 

On the other hand, if the parameter a = 0, then from Eq. (|4.2j) we get the well-known 
formula for a fine structure of the hydrogen atom energy spectrum, for which there is the 
exact solution of the Dirac equation. 

In the absence of the Coulomb interaction (e 2 = 0) the bound states exist for a < 
with the energy levels 



E 



mc 



1 - 



An* 



where the quantum number I*, contained in n*, is determined by equation (|3.4|) for e 2 = 0. 
For the ground state, for n* = + (a/A) 2 — a 2 , one obtains from Eq. ()4.2j) 



E 



mc 



1 + (a/A) 2 



A 



a + y/l + (a/A) 2 - a 2 
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For a = the ground state energy (a < 0) 



E = 



mc 



This quantity plays the role of the particle rest energy, thus, the mean value of the effective 



mass, i.e., 



m 



m 



is smaller than the mass m, which corresponds to the negative parameter a. If parameter 
\a\ ^> A, then the ground state energy E = Hc/\a\ is of the nature of the Casimir energy 
concentrated in the volume ~ |a| 3 . 

Let us now consider the non-relativistic limit c — > oo. We assume that the natural 
length scale for the quantity a is the classical radius of the electron, which, in our problem, 
is the upper limit for a; hence we put a = ae 2 /mc 2 , a < 1. We believe that the numeric 
value of a does not depend on any fundamental constants, and that, ultimatelly, this value 
is the initial characteristic of the particle. Now the formula for the energy takes the form 



E = 



mc 



1 + (a/n*) 2 



£)VJi + £)V*> 




We expand this expression in power series over a up to a 2 : 

-,4 



E = mc 2 



me 
2h 2 n 2 



me 



2h 2 n A 



[l-af 

q2{1 ~ ^ [j+T/2 {1 + 4 (1 + S/3) 



(4.3) 



(4.4) 



where n = n r + I + 1 is the principal quantum number. For a = we have the familiar 
fine structure formula. As we can see, the dependence of the mass upon the coordinate 
also deforms the non-relativistic term (i.e., the Bohr formula) in the way as if the particles 
mass had been substituted with m(l — a) 2 . For a > this effect has the same sign as the 
correction accounting for the finiteness of the nucleus mass. The degeneracy of terms, in 
particular of S1/2 and P1/2, holds because we have the specific dependence of the mass on 
r. If, in expansion (1.1), we also left the next terms, this degeneracy would be removed. 



5 Conclusion 

We have found the exact solution of the Dirac equation for a particle with the position- 
dependent mass, which might be usefull in the study of the corresponding non-relativistic 
problem as a reference result. The next terms, which we have neglected in this work (in 
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particular the dipolar one) and which are responsible for the super-fine structure of the 
energy spectrum, can be taken into account by means of the standard perturbation theory. 

Let us return to the ordering problem in the kinetic energy operator in the Schrodinger 
equation when the mass of the particle depends on the coordinate. We have already 
seen that the corrections from this dependence for the non-relativistic limit are due to 
the fact that in the mass this dependence occurs to be ~ 1/c 2 . The mechanism of 
the apperance of the abovesaid correction is quite simple. The particle rest energy is 
m*c 2 = mc 2 (l + ae 2 /mc 2 r) = mc 2 + de 2 /r. In the linear approximation, according to per- 
turbations theory, this energy is (mc 2 + d(e 2 /r)) = mc 2 + de 2 /a^n 2 = mc 2 + ame 4 /fo 2 n 2 , 
where = h 2 /me 2 is the Bohr radius. This expression, together with the zeroth approxi- 
mation for E, in accordance with the Bohr formula, yields mc 2 — me 4 (l-2d)/2h 2 n 2 , which, 
in the linear approximation in a coincides with formula (|4.4|) . Otherwise said, only on 
condition that the dependence of the mass upon the coordinate is ~ 1/c 2 , we deform the 
non-relativistic expression for the energy levels. At the same time, it is irrelevant how we 
order the momentum operator with the inverse mass in the kinetic energy operator because 
the corrections are of the order of 1/c 2 . 

Thus, the issue of the form of the kinetic energy in the Schrodinger equation for the 
position-dependent mass has, in the general case, quite a limited sense. If this dependence 
arises in the non-relativistic result of the reduction of the many-particle problem 

to a single-particle one, or as problems in the curved space, the specific ordering of the 
operators appears rather naturally in each problem. As a consequence, the energy levels, 
which depend on some ordering parameter, differ from problem to problem. This has been 
shown in Ref. where the Schrodinger equation with the Coulomb potential is solved 
with some other dependence of the mass on the radial variable. Yet it is important that 
for highly excited states, i.e., for large values of quantum numbers, the dependence of 
the energy on the ordering "parameter" is, probably, insignificant, and the energy can be 
found using the quasi-classical Bohr-Sommerfeld quantization method. For instance, for the 
dependence of m* on r given by Eq. (|2.2|) the direct solution of the Schrodinger equation, in 
which \j\fmF stands leftwards and rightwards of the momentum square, obviously yields 
the same energy levels as the Bohr-Sommerfeld method. Finally, there is also another 
argument in support of our theory, which follows from the results of Ref. ^3] , that for large 
quantum numbers the energy does not depend on the inverse mass and momentum ordering 
"parameter" . 

The author is very grateful to Volodymyr Tkachuk for insightful discussions on this 
problem. 



8 



References 

[1] Geller M R and Kohn W 1993 Phys. Rev. Lett. 70 3103 
[2] Serra L and Lipparini E 1997 Europhys. Lett. 40 667 

[3] Barranco M, Pi M, Gatica S M, Hernandez E S and Navarro J 1997 Phys. Rev. B 56 
8997 

[4] Arias de Saavedra F, Boronat J, Pollas A and Fabrocini A 1994 Phys. Rev. B 50 4248 

[5] Bastard G 1988 Wave Mechanics Applied to Semiconductor Heterostructures (Les Ulis: 
Editions de Physique) 

[6] von Roos O 1983 Phys. Rev. B 27 7547 

[7] Levy-Leblond J M 1995 Phys. Rev. A 52 1845 

[8] Cavalcante F S A, Costa Filho R N, Ribeiro Filho J, de Almeida CAS and Freire V 
N 1997 Phys. Rev. B 55 1326 

[9] de Souza Dutra A, Hott M and Almeida CAM 2003 Europhys. Lett. 62 8 

[10] de Souza Dutra A, Almeida CAM 2000 Phys. Lett. 275 25 

[11] Dekar L, Chetouani L, Hammann T F 1999 Phys. Rev. A 59 107 

[12] Green H S 1965 Matrix mechanics (P. Noordhoff LTD-Groningen, the Netherlands) 

[13] Vakarchuk I 2004 Quantum mechanics, 2nd ed. [in Ukrainian] (Lviv University Press) 

[14] Quesne C, Tkachuk V M 2004 J. Phys. A 37 4267 



9 



